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Abstract 

The transverse momentum distribution of prompt photons coming from the 
very early phase of ultrarelativistic heavy ion collisions for the RHIC and LHC en- 
ergies is calculated by means of perturbative QCD. We calculate the single photon 
cross section (A + B —* 7 + X) by taking into account the partonic sub processes 
q + q — * 7 + 5 and q + g^j + qas well as the Bremsstrahlung corrections to 
those processes. We choose a lower momentum cut-off ko = 2 GeV separating the 
soft physics from perturbative QCD. We compare the results for those primary 
collisions with the photons produced in reactions of the thermalized secondary 
particles, which are calculated within scaling hydrodynamics. The QCD processes 
are taken in leading order. Nuclear shadowing corrections, which alter the in- 
volved nuclear structure functions are explicitly taken into account and compared 
to unshadowed results. Employing the GRV parton distribution parametrizations 
we find that at RHIC prompt QCD-photons dominate over the thermal radiation 
down to transverse momenta kx ~ 2 GeV. At LHC, however, thermal radiation 
from the QGP dominates for photon transverse momenta kx < 5 GeV, if nuclear 
shadowing effects on prompt photon production are taken into account. 



Frankfurt, February 1, 2008 
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1. Introduction 

In recent years a lot of effort has been made, on the experimental as 
well as on the theoretical side, to investigate the physics of the quark-gluon 
plasma (QGP) fl]. By doing so, one hopes to gain insight into the state of 
hot and dense matter created in a heavy ion collision and, in a more general 
manner, to learn about the evolution of the early universe which is believed to 
having passed this state shortly after the big bang. During the time when the 
matter, produced in the collision of two heavy nuclei at very high energies, is 
in the quark-gluon phase, particles stemming from the interaction between 
the plasma constituents will arise. By detecting the produced particles one 
hopes to gain knowledge of the physics of the QGP. Those signatures for the 
plasma can be dilepton production, J / ip suppression and photon production 
(for a review, see pf). It is important to know the contribution from col- 
lisions between the initial nucleons, i.e. from a non-plasma source, giving a 
background contribution to the thermal yield. 

In this publication we consider prompt photon production via the QCD pro- 
cesses q + q — > 7 + g and q + g — > 7 + q and by their Bremsstrahlung 
contributions up to 0(af) at midrapidity to investigate by which amount 
nuclear effects as shadowing contribute to the hard process. We take the 
formulas in leading-order and simulate the higher order terms by a K-factor. 
In order to get information about the primary 7's one needs reliable infor- 
mation about the infrared dominated quantities entering the cross sections, 
i.e. the parton distribution functions. The uncertainties appearing in con- 
nection with the QCD cross sections in nuclei are twofold. Not only one has 
to deal with the low- a; behavior of the parton distribution functions but also 
it is unavoidable to account for modifications of those distributions in nuclei, 
such as shadowing corrections, especially when considering heavy nuclei such 
as Pb and Au. 

We consider photon transverse momenta in the range hr = 2... 5 GeV which 
is relevant for the QGP (it was suggested || that in the region kx = 2... 4 
GeV the thermal photons from the QGP may dominate whereas the hard 
QCD 7's should dominate the larger kx region). 

The hard collisions take place at some time r ~ l/&o ~ 0.1 fm/c. The typi- 
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cal time scale for the QGP formation requires knowledge of the dynamics of 
the formation and equilibration of the plasma. We therefore employ two sets 
of initial conditions for the temperature and the proper time in the hydro 
calculations to obtain a lower and an upper bound for the thermal yield. 



2. Inclusive single photon production in QCD 

We first consider the cross section of single 7 production on the nucleon nu- 
cleon level. The general form of a factorized inclusive cross section with large 
momentum transfer kx can be written as M 

a(k T ) = H ^f 2 ®f 2 + (l/k T )H 1 ®f 2 ®f 3 + (l/k 2 T )H 2 ®f^f 4 +O(l/k 3 T ) (1) 

where the process at the nucleon nucleon level is 

N(P, S) + N'(P', S) -> 7 (A;) + X (2) 

as illustrated in figure [I]. The convolution is formulated in terms of the mo- 
mentum fractions x and y. The f n are infrared dominated nonperturbative 
matrix elements of twist n and the Hi are the perturbatively calculable co- 
efficient functions. In our approach we will work in LO with a K-factor for 
the higher orders in a s and at the twist- 2 level. The next twist to contribute 
in an unpolarized process would be twist-4 which is suppressed by 0(1/ k^) 
compared to leading twist and is therefore neglected. Matrix elements of 
twist-3 contribute only to polarized processes, such as single-spin asymme- 
tries u §. 

On the twist-2 level the cross section is 

E C ^-{AB - CX) = ± £ / dx a dx b f a A {x a ) f b B {x b ) (3) 

s — -(ab — > cd) 8(s + t + u) 
dt 

where the sum runs over all partonic subprocess including quark- ant iquark 
annihilation into two photons and into photon and gluon and the QCD Comp- 
ton process. As stated above, we are here only interested in the single photon 
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processes. Reactions involving two final state photons, such as q + q — » 7 + 7, 
are neglected since they are suppressed due to the additional electromagnetic 
coupling constant. Some comments on Bremsstrahlung will appear later be- 
cause one gets a large contribution from those processes. 
Since the above formula applys for the twist-2 level, the Ja,b are the well- 
known parton distribution functions of definite twist. For our calculations 
we used the newest version of the Gliick, Reya and Vogt parametrisations 
|J. The main reason for chosing the newest version is the modifications at 
small x measured at HERA ]7| which show a steep rise in the proton struc- 
ture function F 2 (x,Q 2 ). The photons at midrapidity typically come from 
momentum fractions x ~ 2ko/y/s. Therefore, at the high energy colliders 
treated here (RHIC and LHC) one clearly is in a region where new data on 
parton distributions at low x are important. Another interesting property of 
the '95 GRV parametrisations is the asymmetric sea taking into account the 
violation of the Gottfried sum rule, which defines the difference between the 
proton and neutron structure functions ||: 

Igsr = f-W-K) 
Jo X 

If 1 2 f 1 - 1 

j dx(u v — d v ) H — / dx(u — d) = — if u = d (4) 
3 Jo 3 Jo 3 

The early analysis (prior to 1992) assumed a flavor independent light-quark 

sea. The NMC collaboration, on the contrary, found || 

£(0.004,0.8) = 
r - 8 dx 

/ — (F 2 P - F") = 0.227 ± 0.007 (stat.) ± 0.014 (syst.) (5) 

J0.004 X 

implying d>u. 

The (unpolarized) cross section for the prompt photon production in a nu- 
cleus nucleus reaction is given (without the exchange term) to lowest order 
by [0 

da AB 

E k — (A + B^ 1 (k)+X) = 

- f 1 dx a f t A (x a )ff(x b )^^^(a + b^c + d) (6) 

7T Jx- X a -X l dt 
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where fi(x a ) and /_• (xb) are the parton distributions in the nucleus. The 



various hard scattering functions can e.g. be found in JIT 

( -t 



do 
di 



■naa s e 2 



s 



u . 



do Traa s e z q 8 
~dl ~ P 9 

for the processes q + g — ► 7 + g and g + g — » 7 + g, respectively. 
The partonic variables at midrapidity are given by 



(7) 
(8) 



Xb 



XgX T 

2x a — Xt 



XT 



2kn 



s, x a = x b = k T /Vs, x 



XT 



2-x T 



(9) 



s = x a XbS,u = —XbXTs/2,i = —x a XTs/2 (10) 

We also explicitely take the 0{a 2 s ) corrections to the annihilation and Comp- 
ton graphs, i.e. Bremsstrahlung contributions, into account. Even though one 
has a 0(a 2 s ) contribution, these processes already contribute to 0(a s ) be- 
cause the fragmentation function enters with a factor ln(Q 2 /A 2 ). 
In that process a photon is radiated off a final state quark, or in other words: 
the quark fragments into a photon carrying a fraction z of the initial mo- 
mentum. The fragmentation function is given in leading-log approximation 
by [0 

(11) 
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The result for the Bremsstrahlung contributions is given by JL1 



Ei 



do 



AB 
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n aa 2 „ , k% 1 f 1 , 4 . . . . 



x T 
dx. 



x 



Vt/(2-Vt) X a 

*l + (Vt/2) 2 



F 2 (x a ,A)(xG(x b ,B) + -Q(x b ,B 



x. 



- y T /2 

+ (x a x b , A <-> B 



(12) 



where x b = x a y T /(2x a — y T ) and Q(x) = xY,f[qf{%) + Qf(x)]. The graphs 
giving rise to the different hard parts are shown in figure [| Notice that an 
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expression (xG(xh,B) + ^Q(xb,B)) appears. This so called effective parton 
distribution is due to the different values of the following subprocesses: 



qq^qq : qg -> qg : gg^gg^l : C A /C F : (C A /C F ) 2 . (13) 

One therefore is able to describe the whole process in terms of a single effec- 
tive subprocess with an effective parton distribution 

Q'(x) = xG(x) + ^Q(x) = xG(x) + ^fx + g/x)] (14) 

with C F /C A = 4/9. 

In the next section we will discuss modifications to the cross section due to 
nuclear effects. 

3. Modifications of the cross section: nuclear shadowing 

For the case of prompt photon production in a nucleus nucleus reaction 
one has to deal with nuclear parton distributions entering the cross section, 
xG A (x, Q 2 ) and F A (x, Q 2 ). It is well known from eA, \iA and pA scattering 



pq , P~3| , |14}1 that over the whole x-range, one finds ff~(x, Q ) ^ Aff (x, Q ) for 



the parton densities. Here one has nuclear shadowing, anti-shadowing, the 



EMC effect and Fermi motion, depending on the increasing x-region [[15 
In the present case we are mainly interested in the small-x region, i.e. in 
shadowing. At midrapidity one typically probes the correlators down to mo- 
mentum fractions x ~ kT\fs/(l — kxy/s)- Thus in the interesting k F range 
at RHIC one reaches values x ~ 0.01 and x ~ 3.4 ■ 10 -4 at LHC. 
The shadowing effect described here is a depletion of the parton densities in 
the nucleus. This effect can be understood in terms of parton-parton fusion 
in the infinite-momentum frame If the longitudinal wavelength of a par- 



ton exceeds the contracted size of a nucleon (or the inter-nucleon distance) 
inside the Lorentz contracted nucleus, partons originating from different nu- 
cleons can fuse, resulting in a depletion of the parton densities at smaller x 
and in an enhancement at larger x. One finds from 1/xP ~ 2R n M n /P that 
shadowing shows up at values x ~ 0.1. A saturation of the shadowing effect 
can be expected when the longitudinal parton wavelength exceeds the size of 
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the nucleus. 

One also expects a shadowing effect from a transverse overlap of the partons 
as is seen in the free proton. Here one finds that, for sufficiently small values 
of x and/or Q 2 , the total transverse area occupied by gluons becomes larger 
than the transverse area of the hadron. This happens when xG(x) > Q 2 R 2 , 
with a transverse parton size 1/Q 2 (Q 2 = —q 2 ) and a proton radius R. Hence, 
one finds a depletion of the gluon and sea distribution at values x < 0.01 as 
shown in figure |3|. 

It is obvious that shadowing depends on the one hand on the mass number 
of the nucleus and on the other hand on the x and Q 2 values. The Q 2 de- 
pendence of the ratio (l/A)ff-(x, Q 2 )/ff(x, Q 2 ) is determined by modified 
DGLAP equations (first suggested in [T7| and later proven in flEfl ) taking the 
fusion processes into account |D| . 

In the nucleus those fusion processes are due to an enhancement of the fu- 
sion of gluon ladders stemming from independent partons, whereas in the 
free nucleon, fusion processes from non-independent partons dominate. 
It was shown [19] that the modifications to the usual DGLAP equations 
only affect the ratio by ~ (6 — 7)%. Because the detected photons are in a 
quite narrow range of fcy where k\ ph Q 2 , for simplicity, a parametrization 
for both G A (x, Q 2 ) and F^(x,Q 2 ) is used. The first Ansatz |Tj| assumes a 
similar size of the shadowing of F^(x, Q 2 ) and of the gluon distribution at 
the initial scale Qo = 2 GeV. The evolution of F 2 (x,Q 2 ) is quite moderate, 
whereas the gluon shadowing vanishes faster with increasing Q 2 . This can 
be understood in terms of the momentum flow into gluons at higher Q 2 as 
predicted by the DGLAP equations. 

The parametrisation of the (l/A)F A (x, Q 2 )/F 2 p (x, Q 2 ) ratio at the input 
scale, Qq = 2 GeV, given in fL9| , is shown in figure £|. Nevertheless, the 
same parametrisation for both distributions is used, so one does not need to 
solve the modified DGLAP equations for the narrow kx range. It is clear 
that thereby the shadowing effect is slightly overstimated when one tries to 
be consistent in one Ansatz. 

Due to the electromagnetical neutrality, the Rg{x, Q 2 )-ratio is not accessible 
with deep inelastic scattering processes. A second Ansatz |19| with a much 



7 



stronger shadowed gluon distribution (Rg{x, Q 2 ) is roughly 45 % smaller than 
Rp 2 at x = 0.001) has been presented at the initial value Q = Q = 2 GeV. 
The neglect of the Q 2 dependence of Ra(x, Q 2 ) is absorbed by the systematic 
uncertainty in the Ansatz. But there definitely remains an overestimate due 
to the neglect of the Q 2 dependence of Rf 2 (x,Q 2 ) which changes by ~ 9% 
as Q varies from 2 GeV to 5 GeV at x — 0.001. Therefore, our overestimate 
alters the result towards larger fey by ~ 10%. 

We show the effect of shadowing for the two energies = 200 AGeV and 
y/s = 5.5 ATeV. The shadowing effect increases as the energy increases which 
can easily be understood in terms of the momentum fraction: as one typi- 
cally probes the correlators (here the parton distribution functions entering 
the hand-bag graph) at values x ~ 2kT/y/s at y = 0, this formula also im- 
mediatly explains the drop of the ratio towards larger transverse momenta. 
The decrease of the shadowing correction with increasing can also be un- 
derstood in terms of the lower integration boundary x^ in in Eq. (^) which 
increases with the transverse momentum of the photon. 
Each cross section is computed on the nucleon-nucleon level as da pp / dy. Then 
this quantity is transformed to the number of events per rapidity in a central 
nucleus nucleus reaction by multiplying with the nuclear overlap function 
Taa(P) at zero impact parameter: 

^ = T^(b = 0)^ . (15) 
dy dy 

One should not forget that the results for central collisions are about four 
times larger than those for collisions averaged over all impact parameters. 
Numerically one finds that T AA (b = 0) « A 2 /ttR a , 

TpbPb = — r, T AuAu = — - . (16) 
mb mb 

This procedure is equivalent to neglecting nuclear shadowing effects. To 
explicitly take into account those effects one has to multiply each parton 
distribution with the parametrized ratio Rf 2 = (1 / 'A) / . 
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4. The hydrodynamical calculations 

For comparison we have also computed the transverse momentum spectrum 
of thermal hard photons. As a model for the time-evolution of the ther- 
malized system we assume a three-dimensional hydrodynamical expansion 
with cylindrical symmetry and longitudinal boost invariance p0[ . For the 
RHIC energy we employ an initial temperature of Tj = 533 MeV and an 
initial time of tj = 0.1236 fm/c, while for LHC we assume Tj = 880 MeV 
and Ti = 0.1 fm/c. The resulting final pion multiplicity in central Au+Au 
reactions (i.e. initial transverse QGP radius = 6.5 fm) is dN n /dy ~ 1460 
at RHIC and dN n /dy ph 5300 at LHC, if the hydrodynamical evolution is 
isentropic and if the QGP and the hadron gas consist of u, d quarks, gluons 
and pions, respectively pT| . According to present knowledge p2| , from the 
viewpoint of maximum thermal radiation these initial conditions (and the as- 
sumption of full thermal and chemical equilibrium) have to be considered as 
the most favorable ones. Therefore, we have also computed the thermal pho- 
ton spectrum for a less rapid thermalization and a lower initial temperature 
(n = 0.5 fm/c, Ti = 300 MeV, dN n /dy « 1050 for RHIC and r % = 0.25 fm/c, 
Ti = 650 MeV, dN^/dy pzt 5350 for LHC). We assume a freeze-out temper- 
ature of Tf = 100 MeV. The photon spectrum above kx = 2 GeV, however, 
does not depend sensitively on this number, except in the case of slow ther- 
malization and low initial temperature (r» = 0.5 fm/c, Tj = 300 MeV). 

The equation of state is that of an ideal gas of massive n, rj, p, and uo 
mesons below Tc = 160 MeV. For T > Tq we assume an ideal QGP (massless, 
noninteracting u, d quarks and gluons) described within the MIT bag-model. 
The bag constant is chosen such that the pressures of the two phases match 
at T = Tc (-B 1 / 4 = 235 MeV) thus leading to a first order phase transition. 

The number of emitted direct photons per infinitesimal space-time volume 
in each of the three phases is parametrized as |2~3" 



F dN1 5aa s E/T ( 2M2E \ 

E is the photon energy in the local rest frame. In our hydro calculations 
we fix a s = g 2 /4:TT = 0.3. A logarithmic temperature dependence of a s does 
not alter the results significantly. Eq. ([T7|) accounts for pion annihilation 
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(7T7T — > /ry), and Compton-like scattering {np — > irj, nr] — > nj) off a p or 77 
meson (in lowest order perturbation theory). In the QGP (2 massless fla- 
vors), alternatively, quark-antiquark annihilation (qq — > #7) and Compton- 
like scattering off a gluon (q,q + g — > <? + 7) are considered. In view of the 
uncertainties of the initial conditions and the time evolution of the tempera- 
ture, which enters the thermal photon production rate exponentially, we have 
not attempted to include higher order corrections, hadronic formfactors, or 
additional processes (e.g. a\ — > 717 decays etc.). Finally, the thermal photons 
produced at each individual space-time point are summed incoherently to 
obtain the spectrum of photons emitted in a nuclear collision. 



5. Results 

We have calculated the transverse momentum spectrum (at midrapidity) of 
prompt photons, d 2 N 1 / fadfady produced in the very early phase (l/k « 0.1 
fm/c) of an ultra-relativistic heavy ion collision via perturbative QCD. We 
take the transverse photon momenta in the range = 2... 5 GeV. The 
lower bound is chosen such that perturbation theory is still valid and one 
does not have to take into account matrix elements of higher twist, such 
as (^P\F +a ijj'y + F^il)\P\ . We use a rather moderate estimate for the higher 
order terms, namely K=1.5 for RHIC and K=1.0 for LHC We compare 
those prompt photons to the thermal radiation, calculated as outlined above 
(cf. figures |7| and |9|). In the relevant range we find that the photons from 
Bremsstrahlung dominate over those coming from the QCD Compton process 
by a factor « 1.6[ln(l/a;y) — 1.1] which itself dominates over the annihila- 



tion process as also found in pll and [21]. In our calculation we find that 



the Bremsstrahlung contribution at RHIC lies above the Compton contri- 
bution in the whole fey-range whereas in [D] one had a cross-over point at 
kx ~ 3.5 GeV. This different behavior should be due to the parton distri- 
butions (Duke/Owens set I) used in [IT|. Also, for the parton distributions 
employed by us, prompt photons dominate over the thermal radiation in the 
entire transverse momentum range at RHIC, cf. figure 0. The inverse slope 
1/T = —d/dk? In {dN / hrdhrdy) of the prompt QCD-photons in the interval 
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k T = 2... 5 GeV is « 500-1000 MeV at RHIC and 500-1100 MeV at LHC, de- 
pending on kx- Hydrodynamics predicts that the inverse slope of the thermal 
photons can hardly exceed 500 MeV at RHIC and 700 MeV at LHC, even 
if very rapid thermalization times n ~ and high initial temperatures 
Tj « 3 (-Bt) are assumed. The effect of the nuclear shadowing correction is 



quite moderate at the RHIC energy as also found in minijet calculations [IS 



Note also (figure that for the second set of initial conditions the hadronic 
phase dominates the thermal photon spectrum, due to the strong radial flow. 
A higher freeze-out temperature would, however, reduce this contribution. 
The effect of nuclear shadowing increases for higher energies (LHC, see fig- 
ure |^) and yields corrections up to a factor of ~ 3.5, and therefore is not 
negligible, especially at smaller momentum fractions. Without shadowing 
corrections one has a cross over point with the thermal yield (r, = 0.25 fm, 
Tj = 650 MeV) at k? « 3.75 GeV. When taking into account the shadow- 
ing corrections one finds the cross over point at larger transverse momenta 
(kx ~ 5 GeV). For the higher initial temperature of Tj = 880 MeV, the cross 
over point is shifted to even larger values of h?- 

Thus, even for the lower initial temperature and larger formation time, pho- 
tons from primary QCD processes in the important range hr = 2... 5 GeV 
are negligible at LHC if shadowing is taken into account. Also, for the initial 
conditions expected at LHC, the thermal spectrum is dominated by photons 
produced in the QGP phase (figure [T~0D . 



6. Summary and conclusions 

In this paper we have computed the number of prompt photons with 2 GeV 
< k T < 5 GeV in a central AA collision at energies of y/s = 200 AGeV and 
yfs = 5.5 ATeV at midrapidity to study nuclear modifications of the naive 
extrapolations from pp results. We investigated different sources for the pho- 
tons and also discussed thermal photon production within a hydrodynamical 
model and various sets of initial conditions. We used the nuclear overlap 
function to derive the multiplicities from the cross sections and assumed cen- 
tral collisions (b=0) which yield multiplicities about four times larger than 
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those resulting from impact parameter averaged calculations. 
With an upper bound for nuclear shadowing and a relatively conservative 
estimate for the higher orders, expressed in the K-factor, at RHIC we obtain 
multiplicities larger than those from the thermalized stage of the reaction. 
For the LHC energy of y/s = 5.5 ATeV we find that the cross over point 
between the thermal photons and the primary photons significantly changes 
towards larger transverse momenta when shadowing corrections are taken 
into account. Even if we slightly overestimate the shadowing effect, one can 
conclude that at LHC the background from QCD processes of the initial 
nucleons is of minor importance in the range kx = 2... 5 GeV but will still 
dominate for the very high transverse photon momenta. 
The inverse slope of the QCD-photons 500-1000 MeV) significantly ex- 
ceeds that of the thermal photons (300-500 MeV) at RHIC and is slightly 
larger for LHC energies. 
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Figure 1: General illustration of the direct single photon reaction. The square 
of Fa and Fb generates the nucleon matrix element and H is the hard part of 
the reaction involving the different partonic subprocesses. 
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(a) 




Figure 2: (a) QCD Compton, annihilation and (b) Bremsstrahlung graphs taken 
into account in our calculations. The wavy lines denote gluons and the dashed 
ones stand for the photons. 
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Figure 3: Depletion of xG(x) in the proton as parametrized by MRS in []24]]. The 
/}_ curve corrsponds to a singular behavior (~ x -05 ) as postulated by BFKL 
analysis whereas the small x behaviors of Sq is taken from Regge anya lysis (~ x°). 
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Figure 4: Ratio (l/A)F^(x, Q 2 )/F.[(x, Q 2 ) at the input scale Q = 2 GeV as 
parametrized in []I9| . 
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Figure 5: QCD Compton photon multiplicity distribution with and without nu- 
clear shadowing for ^J~s = 200 AGeV. 
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Figure 6: Same as figure |] but for y/s = 5.5 ATeV. 
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Figure 7: Transverse momentum distribution of prompt and thermal photons 
produced in a central Au+Au collision at y/s = 200 AGeV. 
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Figure 8: Transverse momentum distribution of thermal photons produced 
the different phases at ^Js = 200 AGeV with different initial conditions. 
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Figure 9: Same as [7] but for ^/s = 5.5 ATeV in Pb+Pb. 
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Figure 10: Same as || but for y/s = 5.5 ATeV. 
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